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FEEDBACK CONTROL LAWS FOR HIGHLY 
MANEUVERABLE AIRCRAFT 

Annual Report for the Period February 1, 1994 to January 

31, 1995 


During this year, we concentrated our efforts on the design of 
controllers for lateral/directional control using |l synthesis. This 
proved to be a more difficult task than we anticipated and we are still 
working on the designs. 

In the lateral-directional control problem, the inputs are pilot lateral 
stick and pedal commands and the outputs are roll rate about the 
velocity vector and side slip angle. The control effectors are ailerons, 
rudder deflection, and directional thrust vectoring vane deflection 
which produces a yawing moment about the body axis. Our math model 
does not contain any provision for thrust vectoring of rolling moment. 
This has resulted in limitations of performance at high angles of attack. 

During 1994-95, accomplished the following tasks for the lateral- 
directional controllers: 

1 . Designed both inner and outer loop dynamic inversion 
controllers. These controllers are implemented using 
accelerometer outputs rather than an a priori model of the 
vehicle aerodynamics. 

2. Used classical techniques to design controllers for the 
system linearized by dynamics inversion. These 
controllers acted to control roll rate and Dutch roll response. 

3 . Implemented the inner loop dynamic inversion and 
classical controllers on the 6 DOF simulation. 

4. Developed a lateral - directional control allocation 

scheme based on minimizing required control effort among the 
ailerons, rudder, and directional thrust vectoring. 
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5. Developed p outer loop controllers combined with classical 
inner loop controllers. 

Dynamic inversion is used to cancel the aerodynamic and inertial cross 
coupling torque about the pitch, roll and yaw axes. Pseudo controls are 
used to generate commands to the actual control effectors which are 
ailerons, rudder and directional thrust vector vane deflection. The 
inputs to the inner loop controller are the roil and yaw rate commands 
and the outputs are the roll and yaw rates so that a 2x2 transfer matrix 
relates the inputs and outputs. Ideally this matrix would be diagonal 
with the non - zero terms integrators. The effects of uncertainties in the 
aerodynamics will result in incomplete cancellation of the system 
dynamics by the dynamic inversion control laws. 

The design of the lateral - directional controllers was considerably more 
complex than that of the longitudinal controllers. One complication was 
introduced by the fact that the lateral handling qualities vary 
substantially more with angle of attack than do longitudinal handling 
qualities. In the case of the longitudinal controller design it was 
possible to find a single short period natural frequency and damping 
factor which satisfied the handling quality specifications for all flight 
conditions. This was not the case for the lateral handling qualities. At 
higher angles of attack a much longer time constant (lower bandwidth) 
is required for level one handling qualities. This resulted in the 
necessity for designing different p controllers for different angles of 
attack. These controllers then had to be scheduled with flight condition. 
Since these controllers differed considerably from one another in terms 
of pole-zero configurations, simple gain scheduling was not useful in 
controller scheduling. After considerable effort, we were unable to 
develop a realistic method for controller scheduling; therefore, we 
developed another approach. 

The complexity and the large differences between the directional 
controllers at various flight conditions was the main obstacle in the 
scheduling. In order to simplify the controllers, first-order inner loop 
classical yaw and roll rate controllers were designed. The bandwidths 
of these controllers were scheduled with flight conditions to meet 
handling quality specifications. Then (I controllers were designed for 
outer loop side slip control. The frequency response of these p side slip 
controllers were similar at various angles of attack and we were able to 
select a baseline p controller to be used at all angles of attack. The 
inputs to this controller were filtered in such a way that at other flight 
conditions, the overall input/output relationships were very close to the 
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p controllers. These filters were simple pole and zero cascades 
scheduled with flight conditions. 


The design of the lateral-directional control laws has taken more time 
than we anticipated; therefore, we have requested an additional years 
extension at no cost extension. Our work during the fourth year will 
include the following topics: 


1 Completion of the design of outer loop control laws using p 
synthesis. 

2. Evaluation of stability and performance robustness of control 
laws 

using real and complex uncertainty models. 

3. Integration of the longitudinal dynamic inversion controllers 
designed previously with the lateral-directional control laws 
currently under design 

4. Evaluation of performance of the integrated lateral-directional 
and longitudinal p/dynamic inversion controllers using nonlinear 
computer simulations. 


In previous studies we did not use the full complement of sensors 
available. In this study we will use all sensors and we will incorporate 
weightings on the sensors so that the control laws will optimize the use 
of the sensors based on their accuracy and frequency response 
characteristics. 

Robustness will be determined by structured singular value analyses in 
which errors are modeled as real or complex depending on their 
physical nature. Variations in aerodynamic coefficients, for example, 
will be represented by real perturbations in the coefficients of the 
equations of motion. Other errors, such as high frequency dynamics and 
sensor errors, will be represented as complex functions as they are 
dynamic in nature and can be modeled as transfer functions 

We have shown in our previous work that linear analyses can be very 
useful but are not adequate to determine important quantities such as 
peak acceleration and side slip during extreme maneuvers. Thus we 
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will test our control laws using nonlinear simulations. Maneuvers such 
as velocity vector rolls and turns at high angles of attack will be 
emphasized. 
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Robust Dynamic Inversion for Control of Highly 
Maneuverable Aircraft 


Jacob Reiner,* Gary J. Balas/ and William L. Garrard* 
University of Minnesota, Minneapolis, Minnesota 55455 


This paper presenls a methodology for the design of flight controllers for aircraft operating over large ranges of 
angle of attack. The methodology is a combination of dynamic inversion andstructured singular value (m> synthesis. 
An inner-loop controller, designed by dynamic inversion, is used to linearize the aircraft dynamics. This inner-loop 
controller lacks guaranteed robustness to uncertainties in the system model and the measurements; therefore, a 
robust, linear outer-loop controller is designed using n synthesis. This controller minimizes the weighted H norm 
of the error between the aircraft response and the specified handling quality model while maximizing robustness 
to model uncertainties and sensor noise. The methodology is applied to the design of a pitch rate command system 
for longitudinal control or a high-performance aircraft. Nonlinear simulations demonstrate that the controller 
satisfies handling quality requirements, provides good tracking of pilot inputs, and exhibits excellent robustness 
over a wide range of angles of attack and Mach number. The linear controller requires no scheduling with flight 
conditions. 


I. Introduction 

T HE objective of this paper is to present a method for design of 
llight controllers that provides desired handling qualities over 
a wide range of flight conditions with minimal scheduling. Accept- 
able stability and performance robustness must be maintained in 
the presence of unmodeled dynamics, uncertainties in the aircraft 
design model, and noisy sensor measurements. 

The aircraft considered in Lhis paper is the NASA high angle-of- 
attack research vehicle (HARV), which is typical of future fighter 
aircraft. It is capable of flight at very high angles of attack and 
has thrust vectoring as well as conventional aerodynamic control 
surfaces. 1 The unaugmented aircraft does not meet handling qual- 
ity requirements and some type of augmentation is necessary. This 
paper considers only the longitudinal control, The controller relates 
pilot longitudinal stick input to the symmetric deflection of the sta- 
bilizer and the longitudinal deflection of the thrust vectoring vanes. 

The control design philosophy is to use an inner-loop, dynamic 
inversion controller and an outer-loop, linear n controller. The dy- 
namic inversion controller linearizes the pitch rate dynamics of the 
aircraft; however, since model uncertainties prevent exact lineariza- 
tion, there will always be errors associated with this controller. A 
simple linear fractional transformation model of these errors is de- 
veloped for use in design of the outer-loop /j. controller. This con- 
troller provides pitch rate following by minimizing the weighted 
//ao-norm of the difference between the actual aircraft pitch rate re- 
sponse to pilot stick inputs and the desired response to these inputs 
as given by a transfer function model based on standard handling 
quality specifications. Thus the outer-loop ^ controller is an im- 
plicit model following design, which provides robustness to errors 
due to the lack of exact cancellation of the pitch rate dynamics by 
the dynamic inversion controller. 

Recently a number of papers have appeared that describe con- 
trollers for a highly maneuverable aircraft. In Refs. 2-5, appli- 
cation of linear multi-input/muUi-output (MIMO) control design 
techniques to this problem were presented. In every case, excellent 
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local performance was achieved, although robustness outside the lin- 
ear region was not guaranteed. For a global design, gain scheduling 
is required. Gain scheduling of MIMO controllers can be compli- 
cated, and an alternative is to use control design methods that directly 
consider the nonlinear nature of the problem. One such alternative 
is dynamic inversion. 

The control of nonlinear systems through the use of their inverse 
dynamics is a topic that has received a great deal of attention in 
recent years. 6-13 Lane and Stengel 11 recommended implementing 
such controllers to improve levels of performance and safety over 
conventional flight controller designs. The idea of applying dynamic 
inversion to highly maneuverable aircraft operating in the post-stall 
regime was presented in Refs. 14-16. The simulations and tests de- 
scribed in the first two references and the Lyapunov analysis in the 
third show that dynamic inversion is a good candidate for control 
of aircraft operating at high angles of attack and /or at high angular 
rates. The main assumptions in Refs. 14-16 are 1) the governing 
equations are known precisely and 2) the aircraft states are mea- 
sured or estimated accurately. If either of these requirements are not 
met, the cancellation of the nonlinear dynamics will not be exact. 
This may have serious consequences since dynamic inversion by 
itself does not guarantee any robustness to modeling uncertainties. 
Moreover, even though the nominal inverted system is designed to 
follow handling quality requirements, it can be shown that when 
the inversion of the aircraft dynamics is not perfect, the closed- 
loop aircraft response may not meet those requirements. 17 Stability 
and performance robustness within the dynamic inversion frame- 
work are addressed by synthesizing a robust, outer-loop controller 
around the dynamic inversion inner loop. The outer-loop controller 
is designed using structured singular value (*z) techniques. 

The proposed technique results in the dynamics from the de- 
sired input to the controlled output being linear and constant over 
[he entire Right regime. Therefore, gain scheduling complexity is 
substantially decreased or even avoided through the use of dynamic 
inversion. The inpul/outpul relations are not identical over ihe flight 
regime since the feedback linearization achieved using dynamic in- 
version is not perfect. However, a single n controller can be syn- 
thesized in conjunction with a linearizing feedback such that the 
closed-loop system is robustly stable and meets a predefined de- 
sired handling quality specification over a significant portion of the 
flight envelope. 

II. Dynamic Inversion 

Dynamic inversion is a technique in which feedback is used to 
linearize the syslem to be controlled and to provide desired dynamic 
response. 18 In this paper, dynamic inversion is used to linearize 
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(he pitch rate dynamics of the aircraft. The equation of motion for 
the pitching dynamics is 


. 2 2\ 

q = j-(r 2 - p 2 ) 


,)rp M a M c 

— + — + r 

•yy 1 yy 


(1) 


where p, q, and r are the roll, pitch, and yaw rates, respectively, 
in body axes coordinates; M a is the aerodynamic pitching moment 
excluding the moment from the aerodynamic control surfaces; and 
M c is the control moment due to aerodynamic surfaces and thrust 
vectoring. Dynamic inversion uses the control moment to cancel the 
dynamics of the aircraft and insert a desired response qj. Ideally, 
this results in a linear input /output relationship from the desired 
input qj to the actual output q. Denoting the measured angular 
rates and estimated aerodynamic and physical data by hats (“), the 
commanded moment is 


M c = ■ 
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Fig. 1 General interconnection structure. 


acceleration equation: 




rp + m a + m^q 


4- (A h, K tvc ) 


If all of the data and measurements are available with no errors, the 
control law defined in Eq. (2) results in q — q,i- 

In the HARV, the commanded moment is implemented by de- 
flection commands to the aerodynamic surfaces <5, and the trust 
vectoring vanes <5 /uc . The allocation of effort between these two 
surfaces can be accomplished by a variety of methods and is the 
subject of recent research. 14-16 The usual procedure is to use aero- 
dynamic control for small maneuvers with thrust vectoring being 
incorporated where needed. In this paper, a simple minimum con- 
trol energy algorithm for allocating effort was used; however, the 
dynamic inversion design procedure will work equally well with 
other control allocation schemes. 14-16 In this paper, the deflections 
are allocated to minimize the cost function 


(.y'-Gfc)! 


(3) 


subject to the constraint 

+ (4) 


where and \U t are the estimated values of the thrust vector 
and stabilizer moment derivatives. Minimization of J results in the 
commanded deflections given next: 




<$„ 
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= K lvc M t 

(5) 

= k,m c 

(6) 


The stabilizer becomes ineffective at higher angles of attack. This 
results in a small value of Ms, , and most of the control torque is sup- 
plied by thrust vectoring. At lower angles of attack, thrust vectoring 
and the stabilizer play more equal roles. The commanded deflections 
are the same as the actual deflections when the actualordynamics 
arc neglected, so that the actual control moment is 

M l = (M 6 , K x + Ms iu , k IVL ) 


- - p t + — j — ~ r t } “ 


(7) 


The term M a is approximated in the dynamic inversion implemen- 
tation using a least-squares representation of the aerodynamic data 


M„ » \pV 2 S rcl C rc 
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m* + m^q (8) 


Inserting Eqs. (7) and (8) into Eq. (1) leads to the following pitch 
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i Ml .2 

-r— {r* - p *) 4- : rp - m a — m^q 4- q t( 
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(9) 


Equation (9) is the actual closed- loop input /output relation af- 
ter applying feedback linearization. From Eqs. (5) and (6), it can 
be observed that if K, uc =■ K, uc and K t — K s , then M if k, + 
M itv( k tvc = 1. Moreover, if all of the estimated data and measure- 
ments are exact, then Eq. (9) yields an exact linear relation q = q d 


III. p. Synthesis 

Feedback linearization described assumes exact knowledge of 
the system. Practically this is not the case. To address robustness to 
measurements and model uncertainties, an outer-loop controller is 
designed based on structured singular value (ja) techniques. 19-23 

The general framework of (j. synthesis, shown in Fig. 1 , is based 
on describing the system as a linear fractional transformation (LFTj. 
Any linear interconnection of inputs, outputs, and commands along 
with perturbations and a controller can be viewed in this context and 
rearranged to match this diagram. The term P represents the system 
interconnection structure, A the uncertainties, and K the control 
law. The vector y is composed of measurement signals provided to 
the controller, u is a vector of inputs from the control law, and i 
and uj are the inputs and outputs of the uncertainly block. Robust 
performance of the closed- loop system is achieved if all of the per- 
formance requirements, given in terms of the weighted H norm, 
are satisfied for all possible plants as defined by the uncertainty 
description A. The A block is normalized to I for the synthesis 
procedure. 

The synthesis problem involves finding a stabilizing controller 
K such that the performance requirements are satisfied under pre- 
scribed uncertainties. The D-K iteration, which is an approximation 
to p. synthesis, is a practical approach to design control systems with 
the goal of robust performance. 23-25 This technique integrates H x 
optimization methods for controller synthesis and the structured sin- 
gular value for analysis. An upper bound of p. is found by scaling the 
closed-loop system with a scaling matrix D obtained from solving 
the p. analysis problem. The problem of robust controller design be- 
comes that of finding a stabilizing controller K and a scaling matrix 
D such that the quantity || DF { (P, K)D~ i \\ ao is minimized. One 
approach to solving this problem is to alternatively minimize the 
preceding expression for either K or D while holding the other con- 
stant. For a fixed £>. an optimal control problem is solved using 
the well-known state-space method. 21 - 26 On the other hand, with 
fixed K , the preceding quantity can be minimized at each frequency 
as a convex optimization in ln(D). 21,23 25 The resulting data of D 
can be fit with an invertible, stable, minimum-phase, real-rational 
transfer function. This D-K iteration process is carried ou» until a 
satisfactory controller is constructed. 


IV. Error Modeling 

We now need to develop a model of the errors implicit in Eq 
(9) to construct our robust outer-loop controller. For purposes of 
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Fig. 2 Linear formulation for dynamic inversion analysis. 



Fig. 3 Magnitude and phase bounds of possible q/g d transfer functions 
at a * 5 (full) and 35 deg (dashed). 


formulating this error model we assume that r = p — r = p = 0. 
The model is further simplified by modeling the input as a single 
uncertain effector. Incorporating these assumptions into Eq. (9) and 
linearizing results in the following simplified closed-loop pitch rate 
dynamics. 

a 

J J 

. fz* 


The set of all possible transfer functions from q to q d is defined by 
varying the stability derivatives M q , M a , and Mj, from their nom- 
inal values. The uncertainty descriptions can include model errors 
in (he measurements as well as uncertainties in the aerodynamic 
coefficients. Since ihe level of uncertainty in the aerodynamic pa- 
rameters was not given, a maximum uncertainly of 20% is estimated 
for each of the stability derivatives. A linear representation of the 
modeling errors is derived via the calculation of a number of transfer 
functions from the input q d to the output q. These transfer functions 
are included in the set of all possible transfer functions that arise 
from the maximum ±20% variations in the stability derivatives. For 
purposes of analysis, each of the three parameters, M,, i = a, q, 
takes on only three values: 0.8 A/,, M t , and 1.2 M,. All possible 
combinations of the linear transfer functions, for the same nomi- 
nal system, are calculated. To cover the operating envelope without 
thrust vectoring, the standard F/A-18 suffers wing rock at angles 
of attack greater than 45 deg and ceases to be a stable flight test 
platform. With the addition of thrust vectoring, the HARV is stable 
and maneuverable up to 70-dcg angle of attack. 1 Linear models at 
four different angles of attack — 5. 20, 35, and 50 deg — arc used in 
the development of an error model. 37 Magnitude and phase plots for 
a — 5 and 35 deg are shown in Fig 3. 

The resulting transfer function model describing the set of pos- 
sible linear input-output relationships between q and q d is defined 
by (1.04 -f 0.21 S 2 )/(s + 1 .2S| ) with -1 <<5| and & 2 < 1. The term 
(5 1 corresponds to pole variations and S 2 corresponds lo gain varia- 
tions in the linearized transfer function models. This set of transfer 
functions can be represented by the following LFT in the control 


Z a (I + Z,)" a] 

M i JUJ 

- M a ct - M ,</) (10) 
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Each z, is connected to w, through 8; for i = 1,2. 


( 11 ) 


V. Initial Outer-Loop Controller Synthesis 

Two outer-loop controllers (based on different performance ob- 
jectives) using i± techniques are synthesized and analyzed. The first, 
denoted as K* t is based on feeding back ihe pitch rate measurement 
only, whereas the second controller, denoted as K 2 , makes use of 
both pitch rate measurement as well as the integral of pitch rate. The 
pilch rate and angle measurements are assumed to be corrupted by 
sensor noise. 

It is desirable to avoid actuator saturation. This can be accom- 
plished by weighting the deflection and deflection rate of the actua- 
tors; however, using the dynamic inversion formulation, the actuator 
dynamics are not directly accessible. Thus a weight on the output 
of the fi controller, the desired pitch acceleration, is used to limit 
the control bandwidth. For purposes of synthesis, the maximum 
deflections are ±10.5 deg (0.183 rad) and ±16 deg (0.279 rad) 
for the stabilizers and TVC, respectively For this flight condition 
\f if - -0.735 1/s 2 and = -1.683 1/s 2 . These conditions 
limit the maximum pitch acceleration to 


\q\ mM =0.735 *0.183+ 1.683 * 0.279 = 0.604 



(12) 


Using this value, the pilch acceleration weight is 

_ 1 1 + (s/bw) 

9 ” 0.604 l + (s/100) [_ s : 


(13) 


where bw is a bandwidth parameter to be selected. 

The pitch rate, handling quality transfer function between the 
stick command and the pilch rate was derived from the idealized 
short period aircraft response. 29 30 This transfer function is 

_ K q [\~(s/ZJ\ (| 

stick (s 1 /^) + (2( t p/a> fp )s + 1 

Recent research on handling qualities at high angles of attack in- 
dicated that pilots prefer to have lower short penod frequency and 
higher damping factors as trim angle of attack increases. Using data 
from Ref. 29, the short period damping factor and natural frequen- 
cies were selected as C J/f = landru,^ = 1 rad/s. The numeraior zero, 
Z w , was selected to be -0.312, the open-loop value at an angle of 
attack of 20 deg. The gain K q was selected to be 0. 1 15 rad/s/in. 
so that a maximum slick defleciion of 2.5 in. would command the 
estimated maximum pilch rate. These handling quality parameters 
provide level 1 response from 0 to 30-deg angle of attack. 29 30 Data 
on handling qualities for higher angles of attack are not yet accessi- 
ble in the open literature. The desired handling quality parameters 
could have been scheduled wilh angle of attack, but this would have 
involved some scheduling of the outer-loop controller. 

The weighting function IV, on the error between the actual pitch 
rate and the pitch rate generated by the handling quality model is 




I +U/1Q0) 
l + (s/3) 


1 

rad/s 


(15) 


The preceding handling quality model and weighting function on 
the pitch rate errors lead to a bandwidth of 4 rad/s. Thus a value 
of bw — 4 rad/s was used in the pitch acceleration weight given in 
Eq. (13) lo limit controller bandwidth. 
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Fig. 4 K l configuration: interconnection structure for robust dynamic 
inversion synthesis in pitch rate. 


In the original flight control system, the pilot stick command was 
prefiltered to avoid unwanted high-frequency slick inputs. This filler 
had the transfer function 


24.65 

stick, - j2 + 6 2gl + 9 86 


(16) 


Since the filter was part of the original system, it was also included 
in the interconnection model along with the weighings and the def- 
inition for the inner-loop model uncertainty given in Eq. (11). The 
interconnection model including sensor noise for the first n con- 
troller synthesis is shown in Fig. 4. 

Relating this interconnection to the standard formulation given in 
Fig. 1 results in the following definitions for the inputs and outputs: 




z - Ui zi <7 4A' 

V — [tV| IL) 2 Stick c (17) 

y = (stick, q m y 

u - q ti 

The resulting controller, denoted K J , was 1 3th order As described in 
the next section, this controller yielded good transcient response, but 
its steady-state response was not satisfactory. A second fi controller 
K 2 , which incorporated the integral of the pitch rate, was designed 
and improved the steady-state response. 

VI. Controller Evaluation and Outer-Loop Redesign 

The performance of the controllers are evaluated using a detailed 
nonlinear simulation of the HARV. 30 This simulation includes the 
six degree-of- freedom nonlinear dynamics, high-order models of 
actuator and sensor dynamics, limits on control surface deflections 
and rates, and aerodynamic coefficients obtained from look-up ta- 
bles as functions of angle of attack and Mach number. The baseline 
lateral-directional controller used is provided within the simulation. 
The flap settings are determined by a Mach number and angle-of- 
atlack schedule. There are no elastic modes of the aircraft included 
in the simulation. 

The pitch rate responses with controller K J implemented are 
shown in Fig. 5a for a doublet slick input. Figure 5a illustrates the 
pitch rale calculated from the nonlinear simulations at two different 
trim angles of attack, a = 5 and 35 deg, and the output of the linear 
handling quality model. Note that all of the short-term pitch rate 
responses are very close to one another; however, as time increases, 
the angle of attack of the aircraft drifts from its steady-stale value 
as shown in Fig. 5b. This is because the steady-state pitch rate is not 
zero. The drift is more pronounced at higher initial angles of attack 
and is due to implementation of a controller that was designed using 
the short period mode only. 

Two linear plant descriptions are used to analyze controller K*. 
The first model was only a short period approximation of the air- 
craft. The second model included the phugoid and the short period 
modes of the aircraft. The estimated parameters used in the dynamic 
inversion arc varied according to the relationships given in Eq. (11 ). 
When analyzing the four-stale model, only the two short period 
stales were fed back for the inversion. The linear models for the 
analyses were at 5, 20, 35, and 50 deg. For each angle of attack, 27 


Fig. 5 Nonlinear simulations starting at a - 5 and 3S deg with 
controller A'*. 


cases were analyzed. This resulted in 108 cases for the short period 
approximation and 108 simulations for the four-state longitudinal 
model. The results of this analysis are as follows: 

1 ) All of the closed-loop models were stable using the short period 
approximation. 

2) With the four-state model, unstable closed-loop modes were 
found. The smallest time constant of instability was about 20 s. 

3) For the short period approximation, a is stabilized by the 
stabilization of q. 

It is well known that integral control improves robustness to pa- 
rameter uncertainties and helps to eliminate steady-state tracking 
errors. The simulation results presented earlier indicated very good 
short-term tracking. The long-term performance was not as impres- 
sive as the short period response, due to the divergence of a. This 
divergence was the result of a nonzero, long-term pitch rate. As it 
will be seen, the addition of an integrator leads to a very robust sys- 
tem with respect to both performance and stability. The new control 
problem formulation, which includes the integral of pitch rate, is 
shown in Fig. 6. 

A pilch rate integration is added to the controller. To enforce zero 
pitch rate error at steady state in the H x framework, a 9 handling 
quality model must be included. This is accomplished by integrating 
the output of the pitch rate handling quality model and weighting the 
difference between the output of this integrator and 9 by W 9 . Syn- 
thesizing an H x controller for ihe formulation shown in Fig. 6 leads 
to an extremely large H x norm due to the addition of the integrator 
on the 9 handling quality model. This problem is solvable, using the 
state-space solution, if the subsystem relating the disturbances w, 
to the measurements y does not lose rank on the imaginary axis. 21 
This means that the subsystem from i v to y should be observable 
and controllable on the imaginary axis. Adding an integrator to the 
handling quality model, without bringing its output to the controller, 
leads to the integrator state being neither observable nor controllable 
at the origin. To overcome this problem, the integrator is approxi- 
mated as 1/(5 +0.0001). The weighting function on the error in 9 
is defined as 


= 28.65 


1 + (5/25) 
I + (5/0.5) 


t 

rad 


(18) 
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Fig. 6 K 1 configuration: interconnection structure for robust dynamic 
inversion synthesis in pitch rate. 


Controller Tianifcr Funcwoii* 



10 1 10» 10 l 10J I0 1 

Frequency | Rid/Sec) 


Fig, 7 Bode plots of controller A*. 


The noise added to 9 measurement is ass umed to be white with zero 
mean and spectral density of 0.0 1 rad/Vrad7s. 

The n synthesis was applied to the modified formulation. Only a 
single D-K iteration step was required to achieve the lowest possi- 
ble fi. The controller K * consisted of 1 6 states. Its transfer functions 
are shown in Fig. 7. The integral nature of the stick to q d transfer 
function is apparent. The robust performance ix that was achieved 
is 1.149 with the robust stability being 0.31 1. The y. plots are 
shown in Fig. 8. The low value of the robust stability n implies very 
good robustness characteristics. This is shown to be true from the 
nonlinear as well as from the linear analyses that were performed 
on the closed-loop system with controller K Since the robust per- 
formance jz value was not less than I , we are only guaranteed to 
achieve 87% of the performance objectives for 87% of the size of 
the uncertainty. 

The same simulations presented in Fig. 5 are performed again with 
controller K* implemented. Figure 9 shows the lime histories of the 
most important longitudinal variables. As in Fig. 5, Fig. 9a shows 
two nonlinear pitch rate responses together with the linear handling 
quality model pitch rale response for the given slick input. The non- 
linear responses compare extremely well with the linear response. 
Comparing Fig. 9a with Fig. 5a shows that the differences in the 
amplitudes between the pitch rates is smaller for the K j controller 
than for the K J controller. The difference is very significant in the 
long-term behavior of the system. Steady-stale errors on the pitch 
rales are eliminated, and the long-term angle-of-atlack responses do 
not diverge. The stabilizer did not saturate in deflection; however, 
the de flection and rale of the TVC vanes saturated for very short 
time periods in some cases. 

The controller based on nonlinear inversion with the /z con- 
troller, which includes the integral of pilch rate, provides very good 



Fig. 8 Upper bounds for robust performance and robust stability n 
for controller K*. 



a) Pitch rates of the 2 nonlinear simulations, handling quality and stick 
Input 



T.m. [Sec| 

b) Angles-of-attock 

Fig. 9 Nonlinear simulations starting at a = 5 and 35 deg with p 
controller K*. 


nominal performance. The following tests were applied to evaluate 
the robustness of the controller; 

1) Addition of a constant to the feedback linearization scaling: 
this test was accomplished by scaling the linearizing feedback by 
+20% or -20%. 

2) Addition of white noise to the linearizing feedback scaling: 
instead of being 1 , the new scaling is ( l + a ) where a is white noise 
with zero mean and standard deviation of 0.2 

3) Six degree-of-freedom maneuvers: a longitudinal stick input 
is applied simultaneously with lateral stick commands that excite 
lateral dynamics. 

All of the tests were performed with noisy measurements 

The second lest is the most severe and is illustrated in Figs. 10 
and 11. This test uses a very high frequency time-varying scaling 
etTor in the linearizing feedback. The white noise model on the 
scaling has unlimited bandwidth. Note the high-frequency noise on 
the actuator rates shown in Fig, 1 1. This is due to the direct teed 
through to the actuators of the lineanzing feedback. This actuator 
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a) Pitch rates of the 2 nonlinear simulations handling quality 



b) Angles-of-atlack 

Fig. 10 Nonlinear simulations starting at a = 5 and 35 deg with the 
relation Torque,. = Torque c (l + noise). 



Fig. 11 Nonlinear simulations starting at a = 5, 20, 35 and SO deg with 
the relation Torque,. =Torque c (l ♦ noise). 


noise has almost no effect on the pitch rate response as seen in 
Fig. 10. 

The third nonlinear lest, not shown in the paper, examined the 
ability to decouple the pilch rate response from nonlinear in roll and 
yaw dynamics. 2 '* A lateral stick doublet of ±2.0 i n. was followed by 
a longitudinal slick doublet of the same amplitude. The lateral slick 
was applied to excite the lateral dynamics. It was observed that even 
at the roll and yaw rales of up to 100 and 25 deg/s, respectively, good 
decoupling was achieved. In all of the decoupling tests the pilch rales 
closely measured the response of the handling quality model and 
almost no longitudinal response was observed when applying only 
the lateral stick. 27 Hence, it can be concluded that the decoupling 
properties of the robust dynamic inversion controller are satisfactory. 

A linear analysis was performed for controller K 2 similar to the 
one described for controller K J. As with controller AT], the short 
period approximations of the aircraft exhibited no right-half plane 
poles. For the full longitudinal approximations, only three unstable 
poles re su lied. The fastest unstable time constant was 400 s com- 
pared with 20 s for controller . This is a significant improvement. 


VII. Conclusions 

The robust dynamic inversion methodology described in this pa- 
per has applicability to a wide range of flight control design prob- 
lems. It provides low-order, robust control laws that require little 
scheduling with flight conditions. The resulting dynamic response 
is excellent over a wide range of angles of attack and Mach numbers. 
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